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We study the dynamics of skyrmions in chiral magnets in the presence of a spin polarized current. The 
motion of skyrmions in the ferromagnetic background excites spin waves and contributes to additional damping. 
At a large current, the spin wave spectrum becomes gapless and skyrmions are created dynamically from the 
ferromagnetic state, resulting in a liquid-like state of skyrmions. At an even higher current, these skyrmions are 
strongly deformed due to the damping and become unstable at a threshold current, leading to a paramagnetic- 
like state. We show how skyrmions can be created by increasing the current in the magnetic spiral state. We then 
construct a dynamic phase diagram for a chiral magnet with a spin polarized current. The instability transitions 
between different states can be observed as experimentally clear signatures in the transport measurements, such 
as jumps and hysteresis. 

PACS numbers: 75.10.Hk, 75.25.-j, 75.30.Kz, 72.25.-b 



Introduction - Skyrmions as topological objects were first 
proposed as a model for baryons 1 1 1 and were first realized 
experimentally in condensed matter systems in quantum-Hall 
ferromagnets. ||2l Skyrmions were also predicted to be sta- 
ble in chiral magnets, where the inversion symmetry is bro- 
ken due to the spin-orbit coupling. |[3]-2] In this context, the 
spins wrap a sphere when one moves from the center of the 
skyrmion to infinity. Hexagonal skyrmion lattices were ob- 
served in the A phase of a metallic ferromagnet MnSi f6^ 
as well as in compounds with a B20 structure |7, 8|. More 
recently skyrmion lattices have been directly imaged in thin 
films of FeosCoo.sSi by Lorentz force microscopy |9|, and 
these lattices were shown to be more stable than those in bulk 
materials. |10, 11 1 The continuing identification of more ma- 
terials with skyrmion lattices suggests that they are ubiquitous 
in magnetic metals without inversion symmetry. 

One of the main themes in the field of spintronics is the 
manipulation of spin texture by electric currents. The main 
player is magnetic domain walls (DWs) and their controlled 
motion by current is under active study. It has been shown 
IHJ that the threshold current to move DWs is extremely high, 
on the order of 10^ - 10^ A/cm^, which is detrimental for 
applications due to the intensive Joule heating. Skyrmions are 
also promising candidates for the application of spintronics as 
information carriers, but their potential use for this purpose 
depends crucially on both controlled creation of skyrmions 
and manipulation of their motion. Skyrmions can be driven by 
a spin polarized current in a manner similar to DWs. Recently 
it was demonstrated experimentally that the threshold current 
to move a skyrmion is 4 to 5 orders of magnitude smaller than 
that required to move DWs lITsmSl . indicating that skyrmions 
could have a tremendous advantage over DWs since the Joule 
heating can be reduced significantly. Thus, an understanding 
of the dynamics of skyrmions is both necessary and timely, 
particularly the stability of the driven skyrmions and how fast 
they can be manipulated. 

In this Letter, we study the motion of a single skyrmion in a 
ferromagnetic (FM) background and the flow of a skyrmion 
lattice in the presence of a spin polarized current. Due to 



the damping, the skyrmion acquires a velocity transverse 
to the current (skyrmion Hall effect). The motion of the 
skyrmion peiturbs the FM background causing the radiation 
of a spin wave. This radiation is enhanced when the veloc- 
ity of skyrmion reaches a threshold and the spin wave spec- 
trum becomes gapless, leading to an instability where multiple 
skyrmions are created and the voltage, which is proportional 
to the number of skyrmions, increases sharply. At a higher 
current, the skyrmion structure is strongly distorted due to 
the damping and finally the skyrmions become unstable, re- 
sulting in a paramagnetic-like state with random and chaotic 
precession of spins. The ratio between the longitudinal and 
transverse voltages jumps at the instability. We also show 
that the skyrmion lattice can be dynamically created from the 
spiral state under current. Based on these results, a dynamic 
phase diagram for chiral magnets with spin polarized current 
is constructed. The onset of the different instability transi- 
tions can be observed from transport measurements in exper- 
iments, including pronounced jumps and hysteresis effects in 
the current- voltage curves. 

Model - We consider a thin film of chiral magnets, which is 
described by the spin Hamlitonian with the Dzyaloshinskii- 
Moriya (DM) interaction in a dimensionless form 1131(51 (161 . 
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where 7ex is the exchange constant, D is the DM interac- 
tion due to spin-orbit coupling ((TTl fTgl, and the last term 
accounts for the Zeeman interaction with the external field 
Hfl. Here n(x,y) is a unit vector describing the direction of 
magnetic moment. The external field is perpendicular to the 
film Hfl - HaZ- The phase diagram of chiral magnets obey- 
ing Eq. ([TJ at temperature T = K is summarized as fol- 
lows. 1 16 1 At //fl = 0, the system favors a spiral configuration 
of magnetic moments. Upon increasing the hexagonal 
skyrmion lattice is stabilized at Ha > 0.2D^/7ex- At higher 
fields Ha > 0.8D-/yex, the feiTomagnet becomes the ground 
state. Both transitions are of the first order. 

In the presence of a spin polarized current, the localized 



2 



moment experiences a spin transfer torque due to the conduc- 
tion electrons. The dynamics of the magnetic moments is de- 
scribed by the Landau-Lifshitz-Gilbert equation i20l - [22l 

d,n - (J ■ V)n - yn X Heff + ad,n x n, (2) 

where the first term on the right hand side (r.h.s.) represents 
the spin-transfer torque with spin current density J and the last 
r.h.s. term is the Gilbert damping with damping coefficient a. 
E3l The effective field is Heft = -d'H/dn = /exV^n - 2DV x 
m-B. We solve [24J Eq. (|2| numerically to obtain the current- 
voltage (I-V) characteristics with the electric field given by 
E = n ■ (Vn X d,n). [25 1 The total number of skyrmions in the 
film is Q - J dr^q(r) with the skyrmion density q — a- {dxT^ x 
dyn)l{An). m, 

Dynamics of a single skyrmion - Let us first consider the 
motion of a single skyrmion. In the absence of damping 
ff = 0, Eq. (|2]i has an exact solution (r - vf) with v = -J, 
where (r) is the structure of the skyrmion in the static case 
J -Q. This solution corresponds to the comoving of skyrmion 
and conduction electrons, therefore the spin-transfer torque 
is zero. Thus without damping, the structure of a moving 
skyrmion is the same as that at rest and its velocity is equal 
to -J. With damping a > 0, the structure of the skyrmion 
deforms. For low velocities, however, we can use the quasi- 
static approximation that the structure of the skyrmion is the 
same as that in the static case. The velocity of the skyrmion is 
given by lEH 

V = -J - arjz X V (3) 

with rj - T]ii - 47r[J dr^{dfjn)^Y^ and ji - x, y. The 
skyrmion acquires a velocity component perpendicular to the 
current due to the damping. The velocity parallel to the cur- 
rent is V|| = -J /{I + a^rf) and the perpendicular velocity is 
Vj_ - JaT]/{l + a^if). The corresponding voltages across the 
skyrmion are V\\ - -Anvy and V±_ - Anv^-. Thus the viscosity 
of the skyrmion is 77^, = J/v - -y/l H- a^rj^, and the Hall angle 
is 0H = tan"'(vj^/v||) = -tan"'(a77). We calculate the Hall an- 
gle and velocity numerically to check Eq. ([3|. The velocity is 
defined as v = s(f) with s the center of mass s - J dr^q(r)r. 
We initialize a skyrmion in the FM background in the film 
and then ramp up the current gradually. At a small current 
(velocity), the numerical results of Hall angle and velocity are 
consistent with the analytical calculations, as shown in Fig. [T] 
The motion of the skyrmion is different from that of a super- 
conducting vortex, where the velocity is perpendicular to the 
current. IIZTll 

When the current (velocity) increases, the skyrmion starts 
to deform, especially when a is large, and thus the quasi-static 
approximation in Eq. ([3| becomes questionable. We calculate 
numerically the I-V curve with a single skyrmion at 7 = as 
an initial state. The resulting I-V curve and dependence of Q 
on J are shown in Fig. |2](a-c). For a small current, the motion 
of the skyrmion is stable, as shown in Fig. |2]^d). The move- 
ment of the skyrmion perturbs the ferromagnetic background 




FIG. 1. (color online) (a) The Hall angle tand^ and (b) the velocity v 
vs a for the motion of a single skyrmion obtained numerically (dots) 
and analytically (lines) at 7 = 0.1. 

and excites spin wave excitations. The effect of the radiation 
of a spin wave is more prominent when the skyrmion is accel- 
erated or decelerated by a time-dependent current or scattered 
by impurities. Let us consider the motion of a skyrmion in the 
FM background, which radiates a spin wave n = z + S with 
(5 <K 1 . The dynamics of 6 is 

i(a) - J ■ k)(5^ = -y (//^ + J^^k^ + iacj/yj 6y + hg^A^J, k), (4) 

i(Lo - J ■ k)6y ^y{Ha + 7exk^ + iaoj/y) 5^ + hsjto, k), (5) 

where /!s_^(w, k) - /!s,^(w, k)6{a> - v ■ k) is the eff'ective field 
induced by the motion of the skyrmion. Maximal radiation 
occurs when the resonance condition Q(k) = v ■ k is satisfied, 
where 

nii.)^j-k± '-f^f\ (6) 

is the spin wave spectrum determined from Eqs. Q and (|5]) 
with hs^ft(a),k) = 0. It requires a threshold current Jr above 
which the resonance condition f2(k) = v ■ k can be satisfied. 
For weak damping a <s 1 where v ^ -J, Jr - y y/HaJai/ (a^ + 
1). The dominant wave vector is kjc - y/ Ha /J ex and ky = 0. 
One typical configuration of is illustrated in Fig. [3] which 
clearly shows the radiation of the spin wave due to the mo- 
tion of the skyrmion. For strong damping, the excited spin 
wave decays quickly because the amplitude of the spin wave 
decreases with the damping constant a. The radiation of spin 
waves thus gives an additional contribution to the damping of 
skyrmion. According to the simulations, this contribution is 
small because the amplitude of the spin wave is weak. The 
radiation of collective excitations due to the motion of driven 
topological objects is well known in other systems, such as 
Josephson junctions (281 and DWs 1291 l30l. 

The spin wave excitation spectrum f2(k) [Eq. (|6|] be- 
comes gapless when the applied current reaches the value 
Jin = 2-yyfHaJ^/(a^ + l). This gapless condition signals 
an instability of the FM state into a helical state whose xy- 
spin component has spiral correlations with wave vector k,n = 
(k^ - y/HJJ^,ky - 0). Simulations show that the am- 
plitude of the spiral xy-component grows until the partially 
saturated FM state is destroyed via proliferation of multiple 
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FIG. 2. (color online) (a) E\\IEj_, (b) electric field E, and (c) total number of skyrmions Q vs the spin current J . The arrows indicate the 
direction of current sweep, (d-f) Configurations of n- in the case of (d) a single skyrmion, (e) multiple skyrmions, and (f) the paramagnetic 
phase at different currents. Here o- = 0. 1. (g) Comparison between analytical (line) and numerical (dots) results for the current J„, at which the 
instability occurs and multiple skyrmions are created, (h) Comparison between spatial structure of a skyrmion at rest and at a large current, 
showing distortion. The vectors in the plots denote the and n,, components, and color denotes the n, component, (i) The velocity v, where 
skyrmions become unstable vs a. 



skyrmions [see Figs. |2](c) and|2](e)]. The presence of a mov- 
ing skyrmion affects the relaxation of the original FM state 
into the skyrmion-rich state in the neighborhood of the mov- 
ing skyrmion because spin fluctuations are locally enhanced 
via radiation of spin waves. The amplitude of the k - k,„ 
spin wave keeps increasing until more skyrmions are created 
around the original moving skyrmion. The new skyrmions 
move in response to the applied current and radiate spin waves 
that create a new generation of skyrmions. This multiplicative 
process repeats, until the system is filled with a finite concen- 




FIG. 3. (color online) Spatial structure of the n, component of the 
spin wave radiated due to the motion of a skyrmion. The skyrmion 
location is labeled and the arrow denotes the direction of skyrmion 
motion. Here a = OA, J = 1.4 and = 0.6. 



tration of skyrmions. The electric field increases in a series 
of steps [Fig. [2](b)] because it is proportional to the number 
of skyrmions Q that increases in the same way [see Fig. |2] 
(c)]. The change of current does not modify the direction of 
motion because the ratio E\\/E^ is independent of J [Fig. [2] 
(a)]. Figure |2](g) shows that our analytical expression for J,„ 
reproduces the cuiTent value that is obtained from our simula- 
tions (i.e. value required to trigger the proliferation of multi- 
ple skyrmions). 

The skyrmion is no longer a circular object under large cur- 
rent and develops a tail opposite to the direction of motion, as 
shown in Fig. |2](h). The distortion is stronger for larger a. 
At an even higher current, the structure of the skyrmion is so 
strongly distorted that the skyrmion becomes unstable. Here 
the relation Ej_/E\\ = -arj breaks down as shown in Fig. |2](a), 
indicating that there is an upper velocity limit for skyrmion 
motion due to the distortion induced by damping. Above this 
limit, the resulting state has random and chaotic precession 
of the magnetic moments and the configuration looks like a 
paramagnetic phase [Fig. |2](f)]. When the current is reduced, 
the system recovers to a skyrmion phase at a cuiTent smaller 
than that at which the skyrmion was lost in the ramping up 
process, indicating the occurrence of hysteresis. As the cur- 
rent is further reduced, the number of skyrmions remains con- 
stant down to y = 0. The velocity v, where the skyrmions 
become unstable decreases with the damping constant a as 
Vi ~ 1 / -\Jl + a^rf-, as shown in Fig. jiji). The distortion medi- 
ated instability of a topological object has been demonstrated 
in other systems, such as Larkin-Ovchinnikov instability of 
the vortex lattice in superconductors [31 1 and Walker break- 
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FIG. 4. (color online) (a) Dependence of the longitudinal electric field £1 on spin current at different fields Ha = 0.2, Ha = 0.6, and Ha = 0.8. 
For Ha = 0.2, the number of skyrmions Q is not conserved at high currents and only an averaged value of Q is given. Arrows denote the 
direction of current sweep, (b-c) Configurations n- for (b) the spiral phase and (c) the skyrmion lattice obtained with Ha = 0.2. The spiral 
phase is not perfect and contains some defects (right part of the plot) after the annealing process, (d) Dynamic phase diagram for J vs H^ of a 
chiral magnet obeying the Hamiltonian in Eq. (|T]l in the presence of a spin current. Here a = 0.1. 



down of DWs ll32l . 

Dynamics Starting from the Ground State Configurations- In 
order to understand how the creation and destruction of the 
current driven states could be observed experimentally, we 
next consider the dynamics starting from the ground state 
configurations. The ground states are prepared carefully by 
numerical annealing from high temperature to zero temper- 
ature, where we find a low field magnetic spiral phase, an 
intermediate field skyrmion lattice, and a high field uniform 
ferromagnetic phase. In the feiTomagnetic phase, the effect 
of the current vanishes since V ■ n = in this case, so that 
E = 0. With fluctuations, the current couples with the fluc- 
tuations of magnetization and creates skyrmions at 7„, where 
the spectrum of spin waves becomes gapless. The I-V curves 
for the spiral and skyrmion phases are shown in Fig. |4|a). 
For a small current, the spiral structure does not move and 
no voltage is induced in the system, as illustrated in Fig. |4] 
(b). This intrinsic pinning of the spiral structure is due to the 
DM interaction. Let us consider a spiral structure along the 
X direction n - (0, cos(qx), sin(^x)), where q is the wave 
number of the spiral structure. The spin-transfer torque is 
7f^(0, -sin(^x), cos(^jic)) = -yn x H, with H, = J^qx/y. 
Thus the effect of the spin transfer torque is equivalent to a 
magnetic field along the x direction. This field tends to align 
the spins along the x direction due to the Gilbert damping 
term, which costs energy as a result of the DM interaction. 
Thus, intrinsic pinning of skyrmions exists in clean systems 
(321, similar to the case of DWs Il34l . 

As the current is increased, the spiral structure becomes un- 
stable and a lattice of skyrmions emerges above a critical cur- 
rent, as shown in Fig. |4](c). At 7 0.34, a small number of 
skyrmions (Q = 5) is created while most of the system retains 
the spiral structure. At J ^ 0.6, the remaining spiral structure 
is converted into the skyrmion lattice. These skyrmions move 
freely under an applied current, which induces transverse and 
longitudinal electric fields. The electric fields are proportional 
to the number of skyrmions, so they decrease with increas- 
ing Ha because the skyrmion density decreases. When the 
current is increased further, the skyrmion becomes unstable 



and the system evolves into the paramagnetic phase discussed 
above. At low fields where the spiral structure is the ground 
state, hysteresis occurs upon decreasing the current, while at 
high fields where the skyrmion lattice is the ground state, no 
hysteresis is observed. From these results, we construct a dy- 
namic phase diagram for chiral magnets in a magnetic field 
under spin polarized current, as shown in Fig. |4](d). There is 
hysteresis when the cuiTent is swept across the phase bound- 
ary, so we show the location of the boundary obtained for in- 
creasing current. The hysteresis indicates that the transition 
between different phases has first order features. 

Discussion - Using typical parameters, we estimate the cur- 
rent density at which the instabilities of skyrmion destruc- 
tion occur. The current in Eq. (|2| is given in units of 
yeM^^/fiB, where ~ lO"* Gauss is the saturation magne- 
tization, y X! 10^ Gauss^^s"' is the gyromagnetic ratio ||21J , 
^ a: 10 nm is the size of the skyrmion ID, e is the elemen- 
tary charge, and /ig is the Bohr magneton. We thus estimate 
the instability occurs at a current of order 10*^ A/m^, which 
is easily experimentally accessible. The fastest possible ve- 
locity for skyrmions is of the order of ^yM, ~ 100 m/s for a 
weak damping a <s; 1 and it decreases as l/^fTTahf. The 
Joule heating could be minimized by using a pulse current. 
To examine the instability and motion of a single skyrmion 
in the FM background it should be possible to create a single 
skyrmion experimentally by applying a local magnetic field or 
a circular current ||35]| . 

We emphasize that the dynamic phase transition at 7 = J„, 
occurs regardless of the existence of a skyrmion in the orig- 
inal (J = 0) FM state. The creation of skyrmions at 7„, can 
also be induced by other relaxation mechanisms, such as spin 
fluctuations induced by coupling to the lattice and conduc- 
tion electrons. Such fluctuations provide the dominant relax- 
ation mechanism far away from the moving skyrmion. The 
dynamical creation of skyrmions by current happens even at 
strong magnetic fields where the ground state is FM at 7 = 0, 
which points to a possibly robust way to create skyrmions ex- 
perimentally. In experiments, this dynamical transition can 
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be observed by transport measurements or neutron scattering. 
Moreover, the resulting skyrmions form a liquid-like state, 
which has not yet been realized experimentally by only apply- 
ing external magnetic fields. This new skyrmion liquid state 
opens a door for the investigation of skyrmions in the disor- 
dered state. 
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